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Abstract 

We introduce a transformation which allows the fermionization of operators of any one- 
dimensional spin-chain. This fermionization procedure is independent of any eventual inte- 
grable structure and is compatible with it. We illustrate this method on various integrable 
and non-integrable chains, and deduce some general results. In particular, we fermionize 
XXC spin-chains and study their symmetries. Fermionic realizations of certain Lie algebras 
and superalgebras appear naturally as symmetries of some models. We also fermionize re- 
cently obtained Hubbard models, and obtain for the first time multispecies analogues of the 
Hubbard model, in their fermionic form. We comment on the conflict between symmetry 
enhancement and integrability of these models. Finally, the fermionic versions of the non 
integrable spin-1 and spin-| Heisenberg chains are obtained. 



May 1998 

LAVAL-PHY-21/98 



cond-mat/9806208 



PACS numbers: 02.90.-Fp, 05.50.-hq, 74.20.-z, TLlO.Fd 

Key words: Fermionization, Spin-chains, Integrability, Hubbard models 



'Email address; zmaassarOphy.ulaval.ca 

tWork supported by NSERC (Canada) and FCAR (Quebec). 



1 Introduction 



Fermionization and its reverse procedure, bosonization, consist of a mapping between a set of 
bosonic variables and one of fermionic variables. Such mappings have long been essential tools of 
theoretical physics. They provide a way of expressing the same theory in two different languages 
and open the door to the use of a large number of techniques to study a particular system. They 
also reveal the interplay between symmetries and boundary conditions, and provide new models, 
in particular when the dimension of space or the underlying lattice is changed. Fermionization 
schemes exist in many dimensions. They take an added importance in two-dimensional conformal 
field theories; see for instance Q. 

For the spin-^ one-dimensional quantum spin-chain a fermionization procedure exists which 
allows the mapping between spin operators and fermionic creation-annihilation operators. This 
Jordan- Wigner transformation |^] has been used in particular to study the polaron model |^] 
and the Hubbard model The fermionic forms appear as models of itinerant electrons used to 
model conduction properties, while the bosonic ones appear as spin-chains. More recently the 
Jordan- Wigner transformation was used in the study of the Bariev model of correlated-hopping 
electrons [|[. However the Jordan- Wigner transformation is limited in scope and has long been 
standing alone as a way to fermionize simple one-dimensional spin-chains. It is one of the aims 
of this paper to provide a general fermionization method. 

We start by recalling the Jordan- Wigner transformation and then give its generalization. 
We show how several variants can be implemented to fermionize operators corresponding to 
any spin-chain, integrable or not. We then review how to fermionize, with the help of the 
Jordan- Wigner mapping, an integrable spin-chain model of considerable theoretical and experi- 
mental interest, the spin-^ Heisenberg spin-chain Hamiltonian. As boundary conditions are not 
preserved by the Jordan- Wigner transformation, we clarify this issue on this simple example. 
The fermionic version possesses non-local boundary terms. It is possible to replace these terms 
with periodic ones without spoiling integrability. This example serves also other purposes. The 
general features of integrability, both in the bosonic and fermionic settings, are introduced and 
quantities and equations which will reappear for other models are defined for this simple model. 
We also obtain some general results by studying this system. We then fermionize three of the 
recently found integrable XXC models and the ■su(4) XXZ model, i.e. the model correspond- 
ing to the i?-matrix of the fundamental representation of su(4). The XXC models are hybrid 
trigonometric ones, with an underlying su{2) structure and some su{n) features. They possess 
symmetries much larger than the trigonometric models built from the fundamental represen- 
tations of su{n). We apply the fermionization at the Hamiltonian level and at the integrable 
structure level, by fermionizing the L- and i?-matrices. While different versions of a given XXC 
model exist, they are known to be equivalent at the original bosonic level. This turns out to 
be false for most fermionized versions. We also show that fermionic realizations of algebras and 
superalgebras appear naturally as the symmetries of the fermionized models. Some symmetries 
are broken only to be replaced by new ones. From the examples studied we infer some general 
results. 

The XXC models at their 'free-point' are known to be the building blocks of multistates 
generalizations of the Hubbard model [^, ^. These generalizations inherit the large symmetries 
of the component models. This is another motivation to study fermionized versions of the XXC 
models. We recall the fermionization of the usual Hubbard model and do the same for two gener- 
alized Hubbard Hamiltonians and their L, ii-matrices. This provides for the first time fermionic 
multiple-species generalizations of the Hubbard models. The Hamiltonian density can then 
be written on a lattice of any dimension and provides possible candidates for superconducting 
models in two dimensions. 



1 



The usual Hubbard model has a C/(l) x U{1) symmetry in both its bosonic and fermionic 
forms. In its original fermionic definition, this symmetry is enhanced to an 5'0(4) group sym- 
metry. We investigate whether a similar enhancement takes place for our generalized Hubbard 
models. 

We then derive the fermionic versions of two models of considerable physical importance, the 
spin-1 and spin-| XXZ Heisenberg chains. We propose for the Hamiltonians slightly modified 
versions which do not contain non-local factors. 

We conclude with general considerations and remarks about the foregoing fermionization 
scheme. Finally we provide a preliminary comparison with a, different, recently introduced 
fermionization scheme for integrable models ||8|. 



2 Fermionization 



The mapping from purely bosonic operators to fermionic ones may seem to require the Hilbert 
space, at a given site, to have as dimension a power of two. This is obvious on general grounds. 
The Fock space of a set of anticommuting Fermi operators has dimension 2'^ where d is the 
number of species of spinless fermions, i.e. the number of pairs of mutually anticommuting 
creation-annihilation operators. However, a local space with arbitrary dimension can be embed- 
ded in the nearest larger space of dimension a power of two. The fermionization is then done in 
the latter space, with an eventual projection on the original space. We shall illustrate this proce- 
dure on various examples. Before going any farther we recall the Jordan- Wigner transformation 
for a chain of spin-^ variables [^. 

2.1 Jordan- Wigner transformation 

Let and be the three Pauli matrices: 






= n n = K n = n _i (1 



The transformation at the site m is non-local and given by 
with 

ttJ = I and 7r^_i = cjf • • • C7^_i (3) 

where I is the identity operator. It is straightforward to check that the new operators are 
fermionic creation-annihilation operators, 

where { , } denotes the anticommutator, {A,B} = AB + B A. One may also have more than 
a single copy of the Pauli matrices at every site. For instance, one may have two copies at 
every site resulting in a local space of dimension 4. This is what happens when one bosonizes 
the one-dimensional Hubbard model. The spin 'up' and 'down' fermions result in two copies of 
the Pauli matrices. The p-copies Bariev model contains p copies of the spin-| representation 
at every site Q. Its fermionization results in p species of spinless fermions. The fermionization 
formulae for an A^-sites chain are given by: 

t _ (1) (1)+ _ (1) (1)- 

t _ (1) (2) (2)+ _ (1) (2) (2)- 

^2m — ■^iV+l'^m-l'^ni , a2m — N+l''^m-l'^m (^5) 

at -TT^i) •••vr^P-^^vr^P^ a - n^^'^ ■ ■ ■ ir^^-^K^^'^ a^^^' 
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where 7t!^_i = a[^^^ ■ • • s = 1, ■.■,p. The vrj^^j^ ensure anticommutation between different 

species: 

{asm, atn} = {4m > An) = ^ > {4m > «tn} = Smn^st (6) 

The Pauli matrices satisfy the su{2) algebra. However, and also satisfy the Clif- 
ford algebra and this is the relevant remark which allowed us to generalize the Jordan- Wigner 
transformation. 

2.2 Generalization 

Let cZ be a positive integer. The dimension 2d Clifford algebra, C^d-, is generated by an even 
number of mutually anticommuting generators 

{7i,7i} = +2'^ii , i,j = l,...,2d (7) 

It is a classical exercise to show that this algebra contains the following 2^*^ linearly independent 
elements: 

I; Ti; 7n7i2 > h<i2; ;7n---7i2d' h < ■ ■ ■ < i2d (8) 

The 2d generators of this algebra can be realized as 2*^ x 2"^ hermitian matrices. Thus all 2^ x 2^ 
matrices can be written in a unique way as linear combinations of products of 7 matrices. For 
d = 1 we recover the Pauli matrices while for d = 2 we get the Dirac 7-matrices. There is a 
special element of C2d which anticommutes with all elements 7*, and squares to one. It is the 
product of the 2d 7-matrices; we take a particular normalization for convenience: 

7 = i""* 71 • • • 72d (9) 

Let ^ 

if = -jin2j-\ ± iiij) , i = 1, d (10) 

We define a generalized Jordan- Wigner transformation by the following mapping between bosonic 
matrices and fermionic creation- annihilation operators of d species of spinless fermions: 

4m = ^m-l7/m ' ^^jm = T^^-llJm ' j = l,:;d (11) 

where 

Tr'^ = l and 7r^_i = 71 • • • 7^-1 (12) 
It is straightforward to verify that the following anticommutation relations hold: 

{ajm, akn} = {a]^, 4n} = ' {^m' = ^jk^mn (13) 

Since 7 squares to one, it is invertible and the above mapping is one-to-one. Let 

rijm = Oj^j^ajm (14) 

be the number density operator for the species j. One can then express the matrix 7 at site m 
in terms of fermions operators 

d 

7m = ]l('^nim-l) (15) 

i=l 

Since njm squares to one 7^ also squares one as needed, and the inverse equations are given by 

m— 1 / d \ m— 1 / d \ 

= n EC^riik - 1) 4™ , 77m = n EC^^lk - 1) %m , i = 1, -.d (16) 
k=l \l=l J k=l \l=l ) 
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Some remarks are in order. One could have chosen to group the 7^ matrices two by two in 
an arbitrary order, with an eventual simple change in the definition of the 7 matrix. However 
there is no loss of generality in the foregoing choice, as no specific choice of representation of 
the algebra C2d has yet been made. 

As we noted earlier, one should not be left with the impression that the Hilbert space at 
every site must have dimension 2'^ to be able to fermionize. Indeed any space can be embedded 
in the nearest higher dimensional space of dimension 2'^. More precisely, the operators at this 
site can be written as linear combinations of E^^ matrices whose only non-vanishing elements is 
a one at row i and column j, with i, j varying from one to the dimension of the space. One first 
embeds the matrices E'^^ in the nearest 2'^-dimensional Hilbert space by adding the required rows 
and columns of zero, and then proceeds with the fermionization. The transformed quantities 
act in the larger Hilbert space but also act in a stable way when restricted to the original space, 
where one has identified which states of the larger Fock space span the smaller one. One obtains 
in this way fermionized quantities for the chain. We shall illustrate this procedure on two simple 
examples. 

The above transformations admit several variants. The spaces at different sites of a chain 
need not be equal or of equal dimension. Indeed, all that is required for all that was said above 
to still hold is the one-to-one mapping at one site and the anticommutation relations. And this 



can be ensured by replacing the string 71 • • • 7m- 1 in equations (11,12) with 

— 11 lm-1 V-L' j 

where 7^- is the special gamma matrix acting on the local space of dimension 2 ' at site i. The 
dimension 2^*' depends on the site. The also correspond to the dimension 2"^*. Thus one can 
fermionize for instance a spin-chain with different representations at every site. 

One can also choose to fermionize selectively, leaving some sites unfermionized while fermion- 
izing others. This is achieved through another change in vr^_i. For the sites where one keeps 
the bosonic variables, the bosonic operators are left unchanged by the mapping: i?j — > For 
the sites fermionized, one uses another version of (^) where the string has been replaced by one 
where the identity replaces the 7-matrix for all the sites not fermionized. One thus has ensured 
that the fermionic operators anticommute among themselves as needed, while they commute 
with the bosonic operators for sites where no fermionization has taken place. In this way one 
gets a chain of mixed bosonic and fermionic local spaces. 

There is also the multiple-copies fermionization, when a local space is a tensor product of 
smaller spaces. The fermionization formulae are just the 7-generalizations of formulae (^). One 
replaces the p copies of Pauli matrices by p copies of gamma matrices. Note that copies of gamma 
matrices at the same site need not correspond to the same dimension, as the space at a given 
site may be the tensor product of spaces of different dimensions. Again, all the above variants 
can still be implemented within the different copies, with copies of gamma matrices at different 
sites possibly having different dimensions. The selective fermionization is still possible. We shall 
give examples of fermionizations with two copies of 7 matrices. Finally, reverse fermionization 
i.e. bosonization, with all its variants, can be implemented by use of the inversion formulae 

(H). 

For the specific calculations to follow we have adopted the definitions of 7-matrices found in 
P]. We now fermionize operators of various spin chains. 
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3 Integrable spin-chains 



We start with a well-studied model, the spin- 2 XXZ spin-chain. Its Hamiltonian as well as its 
integrable structure have already been fermionized. This simple example is a warm-up exercise 
which allows us to introduce the general framework of integrability and general aspects of the 
fermionization procedure. We shall then fermionize the simplest XXC models. Unless otherwise 
indicated, roman symbols correspond to the original bosonic quantities, while calligraphic ones 
correspond to the fermionized ones. 

3.1 The spin-i XXZ chain 

We first consider the Hamiltonian and the (L, i?)-system, and then boundary and symmetry 
issues. 

3.1.1 Hamiltonian and L-matrix 

The Hamiltonian is given by 

m 

= IY. ± <^l+i + 2A a^cT^+i) for x = ±l 

m 

Here we have introduced a twist parameter x. In the notation of [jl^ one is considering the 
model (ni = l,n2 = 1). The parameter x arises from threading the ring-chain by a flux N(j). 
One then have x = e*^'?^(^™') where 

(Pixm) = I A{x)dx (19) 

is the Peierls phase. 

The Jordan- Wigner transformation yields 

^2 = - ^ aj^+ittm + x^^aj„a.m+i - A (2nm - l){2nm+i - 1)) (20) 

m 

where A = —2^^, and Um = al^CLm is the number density operator. Let us also introduce 
another number density operator, 

as this quantity appears naturally farther on. Note that the transformation am < — > ajn inter- 
changes rim and fim- This is the particle- hole symmetry which reflects the fact that the roles of 
the creation and annihilation operators can be reversed. Since operators are independent of the 
choice of a fermionic vacuum a and play symmetrical roles, and it is only natural that as 
well as Um should appear. 

The bosonic and fermionic models, both with their respective periodic boundary conditions, 
are known to be integrable. The fermionization of the L-R system was done in [^. The fermion- 
ization of the integrable structure consists in applying a 'gauge transformation' on the Lax 
matrix at every site. The i?-matrix intertwining two such matrices is then obtained by plug- 
ging the fermionic L-matrix into the RLL equation and by cancelling out the transformation 
matrices. 
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We adopt the notations of [|T0| and consider a transformation for L which differs shghtly from 
the one found in 1^, the difference is however not essential. Let a = sin(7 — A), b = sin A, 
c = sin 7, and -E"^ be the matrix whose sole non- vanishing entry is a one at row a and column 
(3. The L = R matrix is given by: 



Define the operators 



Vm = diag(w„, vj)= [ ^"^ I , Vm = ex.p i ^-^ J2 (^i " 1) 1 (23) 



v„ 

The fermionic £-matrix is obtained as: 

■Cm = Kn+l-^m^^ (24) 

and one easily derives: 

/ arim - ibxfim -icam \ 

^m= [ t --1^-1 (25) 

y Ca^rn laUjn + bx J 

The RLL equation is given by 

R{Xi - A2) L(Ai) L(A2) = L(A2) ® L(Ai) i?(Ai - A2) (26) 
and R satisfies the Yang-Baxter equation 

-Rl2(Al — A2)i?23(Al)-Rl2(A2) = -^23 (A2)-Rl2 (Al )-R23 ( Al — A2) (27) 

Let P be the permutation operator on the tensor product of two spaces: P{x®y) = y0x, where 
X and y are two vectors. Define -R(A) = PR{\), then an equivalent form of equation (|2^ ) is: 

R{Xi - A2) L (Ai) i (A2) =L (A2) i (Ai) i?(Ai - A2) (28) 

1 2 
where L (Ai) = L{Xi) I and L (A2) = I i(A2). The Yang-Baxter equation for R is given by 

iil2(Al - A2)i?13(Al)i223(A2) = i?23(A2)iil3(Al)i?12(Al - A2) (29) 



To find the fermionic version of ( Pq ) one inverts ( [2^ ) to obtain L{Xi) and L(A2), and replaces 
them into equation (|2^). Upon commuting the V^s through, one finds that they cancel each 
other out after leaving phases which modify the matrix R. Equation (p6|) becomes 



7^(Al - A2) £(Ai)®£(A2) = £(A2)®/:(Ai) 7^(Al - A2) (30) 
where the tensor product of two operators A and B is graded according to 

(^®i?),,-fci = (-l)(^»+^W)^(^-)^,fci?,, , P(i)E Z2 (31) 

and 

TZ{X)=F-^R{X)F , F = diag(l,l,i,i) = (diag(l,i))0l2 (32) 
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For ( pSD one has: -P(l) = and P{2) = 1. Note that F is a diagonal matrix defined up to an 
overall normalization. We shall choose throughout this paper Fn = 1. This gives the above F, 
and = I4. 

One can then define a graded permutation operator V whose action on two vectors of well- 
defined parity is given by: V{x®y) = {—l)^^^'^^^y^y®x. In components one has: Vij^ki = 
(^-l)P{i)P{i)SiiSjk. Let 7^ = Vn, and equation (|2|) becomes 

n{\i - A2) C (Ai) I (A2) =1 (A2) C (Ai) 7^(Al - A2) (33) 

1 s 2 s 

where C (Ai) = £(Ai)0 I and £ (A2) = I (E>/3(A2). The matrices 7?. and IZ satisfy the Yang- 
Baxter equations ( [27| ) and (29) respectively. However, while (^) written in components for 7^ 



is unchanged, equation (|29D for 7^ has grading signs due to the graded permutation operator. 

We shall obtain similar results for the fermionization of all the integrable systems considered 
here. The above equations will still hold, with the corresponding matrices R, L, F, and the 
appropriate grading. 

3.1.2 Boundary terms and integrability 

It is well-known that the Jordan- Wigner transformation does not conserve periodic boundary 
conditions. One finds 

'^iv^r = -(2'^2 - 1) • • • (2nAr_i - l)ajvai , a]^af = -{2n2 - I) ■ ■ ■ {2nN-i - l)alaN (34) 
Thus the fermionic Hamiltonian ( pO|) with the JW-twisted fermionic boundary conditions 



QN+i = +(2712 - 1) • • • {2nN-i - l)ai , ajy, ^ = +{2n2 -!)••• (2nAr_i - l)a\ (35) 



is the Hamiltonian equivalent to the bosonic one ( [T^ ) with periodic boundary conditions. If in- 
stead one takes periodic fermionic boundary conditions for (|20|), with ajsf+i = «i and a^j^_^-^ = a\, 
the resulting Hamiltonian is not anymore equivalent to the bosonic one with periodic boundary 
conditions. It is however still integrable. This follows from the fermionic relations (^0|). 

To this end we recall the general tenets of integrability. The transfer matrix, t(A), is the 
generating functional of the infinite set of conserved quantities. Its construction in the framework 



of the Quantum Inverse Scattering Method (QISM) is well known [|l2| , [1^ , 14|. Given an {L,R) 



pair, the trace over the auxiliary space of the monodromy matrix T(A) yields the transfer matrix: 

r(A) = TVo [T(A)] =Tro [Mq Lo^(A) • • • Loi(A)] (36) 

where is the number of sites on the chain and is the auxiliary space. The introduction of the 
numerical matrix M corresponds to integrable periodic M-twisted boundary conditions, with 
M = I corresponding to periodic conditions. This twisting of the periodic boundary conditions 
is different from the one arising from the Jordan- Wigner transformation. For matrices M such 
that [M M, R{X)] = equation (|2|) imphes the RTT relations: 

R{Xi - A2) r(Ai) r(A2) = r(A2) r(Ai) r{Xi - A2) (37) 

Taking the trace over the auxiliary spaces one obtains [r(Ai), r(A2)] = 0. A set of local conserved 
quantities is given by 

/dPlnr(A)\ , , 

Hp+i = { ' ] , P>0 (38) 



A=0 



The Hamiltonians -ffp+i mutually commute and the system is said to be integrable. 
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In the fermionic setting the monodromy matrix is given by 



r(A) = Mo£oiv(A)---/:oi(A) (39) 

This is so because because the Lax matrix is a homogeneous even-parity matrix. The general- 
ization of the QISM framework to the graded context was initiated in ||l^, [l^. Equation (|30|) 



replaces equation (|26D, with the appropriate grading, and equation ( pT]) is replaced by 

7^(Al - A2) r(Ai)®r(A2) = r(A2)®T(Ai) n{\i - As) (40) 

It can be shown in all generality that periodic fermionic boundary conditions correspond to a 
matrix M which implements a supertrace: 

M,, = 5,,(-l)^« , Tr(iVM) = Str(A) = ^(-1)^«^,, (41) 



This shows that Hamiltonian ( pO| ) with fermionic periodic boundary conditions is also integrable. 

We restrict ourselves to M = I for the bosonic case, and to a matrix implementing the 
supertrace in the fermionic setting. The latter matrix corresponds to the grading arising in 
equation (^). For the latter case we omit M and replace the trace by a supertrace in (|3^). 
From now on, unless otherwise indicated, 'bosonic system' means we consider bosonic periodic 
boundary conditions, while 'fermionic system' means we consider periodic fermionic boundary 
conditions. Thus we are considering the bosonic and fermionic systems with their 'natural' 
boundary conditions; as we have seen, these systems are inequivalent. 

An immediate and general consequence of all this is the existence of at least two integrable 
boundary conditions for any integrable system: Periodic and possibly many per iodic- twisted 
conditions through fermionization or bosonization. This applies to the general fermionization 
scheme we introduced. 



3.1.3 Symmetries 

Let us now consider the symmetries of the bosonic spin-^ XXZ model. For arbitrary A the 
operator = J2iLi commutes with the transfer matrix and therefore with all the conserved 
quantities. At the rational point where A = |, the full symmetry is su(2), with the two 
remaining generators being = Y^f=i(yf and a~ = Y^f=i(y^ ■ Their fermionic counterparts 
are given by: 




N i-l \ N 



a- =Y.\\{{2n,-l)\a,, = ^(2n, - 1) (42) 



To show that commutes with the fermionic transfer matrix, we showed and used 

[n^,C^{\)\ = -\H,Cm{\)] (43) 

At A = ^, are cj^ still symmetries of the periodic fermionic system? We have verified that the 
fermionic Hamiltonian does not commute with cj^. Thus, although the rational limit of C and 
TZ is well-defined, the symmetry enhancement at the rational point does not take place for the 
fermionic system. This is in fact a general (negative) result for all models fermionized, whenever 
the rational limit provides an enlarged symmetry for the bosonic system. 

This simple fermionization already illustrate features that will also arise for most if not all 
integrable systems. 
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3.2 The (1,2)-XXC model 

This models naturally appears when one takes the infinite coupling limit of the Hubbard model; 
see for instance |17]. The dimension of the the Hilbert space at one site is three. We therefore 
utilize the C4 Clifford algebra to fermionize this model. 
The Hamiltonian is given by 



(44) 



We took ^ = {1} and B = {2, 3} in the notation of |jT^. In the bosonic setting the two other 
choices are unitarily equivalent to this one [18|. 

There are four possible equivalent embeddings of the 3x3 matrices into 4x4 matrices: 



(1,2,3) 
(1,2,3) 



(1,2,3,4) 
(1,3,4,2) 



(1,2,3) 
(1,2,3) 



(1,2,4,3) , 
(2,3,4,1) . 



(45) 



For instance, the second choice corresponds to letting 

f an 

0.21 






012 
022 


032 



while the fourth choice corresponds to 

/ 











ail 
021 
031 





ai2 
a22 
032 



ai3 \ 

023 



033 J 





ai3 

023 

033 / 



(46) 



(47) 



In the specific representation we adopted for the 7-matrices, the canonical basis (61,62,63,64) 
is sent to the four fermionic states as follows: 



61 



63 








1 



44 1 0) 



62 



4lo) 



64 



/o\ 
1 



vo/ 
^o\ 




V 1 / 



|0) 



(48) 



4|0) 



We choose, for obvious reasons but again without loss of generality, the fourth embedding which 
means that the Hamiltonian and the higher conserved quantities. Tip, automatically annihilate 
all states of the chain which have a double occupancy on at least one site. The third choice 
corresponds to the annihilation of all states which have a vacuum on at least one site. 
The fermionized Hamiltonian density is given by 



H 



mm+l 



2^23 «lm^l»Tt+l + 2:23 aim+iaim)n2mn2m+l 



(49) 



+ {X2i 4m02m+l + 2:24 '3^2m+l'3^2 
- ^COS7 (1 +CmCm+l) 
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where C = (2ni — l)(2n2 — 1) — nin2 is the 'conjugation operator' It is obvious that the 
action of 7i2 is stable on the subspace of dimension 3^. Define the projection operator which 
annihilates any state with double occupancy on at least one site: 

N 

P= n ~ nimn2m) (50) 
m=l 

Ti.2 can be interpreted as acting on the reduced Hilbert space or as pTi.2P acting on the full 
4^-dimensional space. 

The Lax matrix is fermionized with the help of the matrix 

m-l 2 



Vrr 



One easily finds 



where Vr, 



ITT 



exp 



E EK^ - 1) 



(51) 



i=i j=i 



/ /:22 

-Caimn2r, 

\ ca2mnim 



r ■ t 



(52) 



-^33 
£44 



Note that (1 
P(3) = P(4) : 



-a nimn2m - i bi^23 nimn2m + X2A nimn2m 
i(inimn2m~ b nim'>T'2m 
i a nim'^2m b X24 f^lm^2m 

f^lm^2m)^m — ■£m(l ^lm^2m) — ^ 

1. We obtain for the matrix F 



The grading is given by: ^(2) = 0, 

(53) 



F = diag(l, 1, 1, —i, —i, —i, —i, —i, —i) = diag(l, —i, —i) ® I3 
The symmetry generators for the conserved quantities are given by 

— Z^m "'Im = Z^m '^2m 

Recall that the off-diagonal generators are symmetries provided the following constraint hold: 



(54) 



X23 = 3^24 1 10 1- We see in the next section how to show similar commutation relations with the 
fermionic transfer matrix. The symmetries are the generators of su{2) x u(l), and reflect a (par- 
tial) symmetry between the two species of spinless fermions. Because of no-double-occupancy, 
X — Zm ''T'imn2m is a trivial symmetry. Note that, here and below, these fermions can be seen 
as one type of fermion with spin ^. We however reserve the terms 'up' and 'down' for the two 
XX copies used in constructing Hubbard models. 

3.3 The (2, 2)-XXC model: first avatar 

We first write this XXC model with A = {1,2} and B = {3,4}. The fermionic Hamiltonian 
density is given by 



n 



(1) 

mm+1 



— a\^aim+i 
+ a\^a2m+i 

+ '^2m+l'^2m 



X2'i ^2m^2m+l — X{ln2mn2m+1 
X2'i^2mn2m+l — Xl4n2mn2m+1 
X24. nimnim+1 — Xi^nimnim+l 
X2lnimnim+1 — Xi^^nimnim+l 



(55) 



icos7(i + c«cSi; 
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where C^^) = (2ni - l)(2n2 - 1). 

The Lax matrix is fermionized with the help of matrix 



We find 



Vrr. 



where 



exp 



^(1) 

'11 



-a aim02m 

'-22 

C Ojlm'^2m 

ca2mnim 



m-l 2 

i=i j=i 



1) 



-ica2mnim 

'(1) 

-33 



-i a at 



% C CLlm'^2m 

—i a a2m,'^lm 
^44 



(56) 



(57) 



'-'11 

'-'22 

A^) 

■'-'33 



a "-lmra2m - « h {Xi3 nimn2m + Xii nimn2m} 
-a nimn2m - i b {^23 nimn2m + X2A nimn2m} 

~1- 



i-44 



X23 ^llm^2m I 

- = ,} 



i a nimn2m + b \^x^^nimn2m 

ianimn2m + b \^X^inimn2m - X24nimn2r. 

The tensor product of equation ( |30| ) is now graded according to 

P(l) = P(2) = , P(3) = P(4) = 1 

For the matrix F we found 

F = diag(l, 1, 1, 1, —1, —1, —1, —1, i, i, i, i, i, i, i, i) = diag(l, — 1, i) ® I4 



(58) 



(59) 



Again one has = Iig. 

Now consider the symmetries of the transfer matrix. It is easy to check that the following 
commutation relations hold provided all the parameter Xap are equal to each other, but not 
necessarily to ±1: 



{4 



(1)1 
m 

(1)1 

Tl J 



^(1)1 



2m^lm.-> } — [-^0 ' _ 



Let 



\P'lm'^2rm ] 

[ni^n2™,/:«] 



^12 

^34 



m J 

11 



m 



[aimCL2m^^m \ — ~[Eq^ , Cm ] 



['^lm^2m) -^m 



(60) 



E 



m "'2m"'lm 
t 

Z^m ^2mflim, 
Em '^lm^2m 
Em "-lm"-2m 



^21 
^43 
^22 
^44 



Em '^lm'3^2m 

t 

Z^m (^lmO-2m 
Em ^lm"'2m 
Em ^lm^2m 



(61) 



Relations 
system: 



imply the following commutation relations with the transfer matrix of the fermionic 



{Xi2,r«(A)} = 



fX34,r(i)(A)] 



11 ^(1) 



(A)] 



[X33,r«(A)] 



{X2i,r«(A)} = 
[X43,rW(A)] 



22 ^(1) 



(A)] 



[X4^r(l)(A)] 









(62) 



We come to the unusual fact that the zero parity operators X^'^ and X'^^ anticommute with the 
zero parity operator t(A) ! However, a simple proof by induction shows that the logarithmic 



11 



derivatives of the transfer matrix are made up of sums of product of an even number of the 
transfer matrix and its derivatives. This shows that the symmetries X^'^ and X^^ still com- 
mute with the conserved quantities Tip, for p > 1. The anticommutation with exp(7Yi) does 
not of course prevent the simultaneous diagonalization of X^'^ or X'^^ with all the other con- 
served quantities, as the spectral spaces are stable under anticommutation and/or commutation 
relations. The symmetry operators clearly are the generators of su{2) x su{2) x with a 

fermionic realization of the generators of this algebra. 



3.4 The (2, 2)-XXC model: second avatar 

We consider here the (2, 2)-XXC model in another realization, where A = {1, 3} and B = {2, 4}. 
As first pointed out in ||l^, this bosonic model is unitarily equivalent to the bosonic model of 
the preceding section. The orthogonal operator linking the second avatar to the first one is 
constructed as follows. Consider the orthogonal matrix 



u = u- 



1 \ 

10 

10 

\ 1 y 



(63) 



The A^-fold tensor product U^^^ = U ^ ■ ■ ■ (^i U transforms the models into each other, modulo 
redefinitions of the parameters Xai3- This can be seen at the Hamiltonian level or generally at 
the transfer matrix level 



f/Wr«(A) C/W=r(2)(A) 



(64) 



At the level of the L-matrices one has: 

Note that C/(^) is an orthogonal and therefore unitary operator. 
The fermionized Hamiltonian density is now given by: 



(2) 

mm+1 



I J 



^12 '^2m'^2m+l 



Im+l'^lm P32 ^201^2™+! — 2^14 ^201^2™+!] 



''lm+ 



Xl2 a2m+l"2m — 3^34 «2m"2m+l 



^cos7(l + C^')cS|i) 



(65) 



(66) 



where C*^^-* = 2ni — 1. The ^-matrix, the induced grading and matrix F are the same as in the 
previous section. We find: 



£(2) 



'-11 



—caima2m 

'-22 

— Caimn2m 
a«2m"-lm 



-t a a2mnim 
-ica[^n2m 

,{2) 
-33 

,t n 
>'lm^2m 



-ica\ 



-icaimn2r. 

^ ^ C'2m^im 

I C a.2m'^l'Ti 
.(2) 



(67) 



44 
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(2) 

C\i = a nimn2m - b {xi2 nimn2m + i Xu nimn2m} 

Cf^ = i a ni^n2m - b {X32 nimn2m - i 2:34 ni„n2m} 
£44-* = ianimn2m + b^x^lnimn2m + i x^lnimn2mj 

This matrix is not fundamentally different from the one in the preceding section. 

Unlike the fermionic version of the first realization, the symmetries which survive fermion- 
ization are only the diagonal ones: the X**, i = 1,...,4. They generate u{l) x x 
There are no constraints on the parameters Xap- The off-diagonal bosonic symmetry generators 
become non-local in the fermionic setting and it is easy to check that they do not commute with 
the fermionic Hamiltonian. This breaks the bosonic su{2) © su{2) © n(l) to the three u{l)'s. 

The obstruction to symmetry comes from the periodic fermionic boundary conditions. Lo- 
cally the two fermionic versions are however equivalent. At the level of the C matrices one 
has: 

Cg^ = (Vm+iUoUmV^l,) {V^UoUmV^') (68) 
The operator Um is non-local: 



m— 1 



l^m = n2m ~^ n ~ l)(2'^2i " 1) j n2m(a{„ + aim) (69) 
At the Hamiltonian level the correspondence is: 

n^2n+l=i^N---Kin^2n+li^l---i^N , m = l,...,iV-l (70) 

However this equivalence breaks down for the last term: TCni- Thus the two fermionic models, 
despite sharing a common bosonic antecedent, are not equivalent. This is due to the periodic 
fermionic boundary conditions. 

This is not the end of the story however. By fine tuning the twist parameters one can recover 
a large superalgebra symmetry. Integrability is preserved since it does not depend on the values 
of these parameters. We take 

X = Xi2 = X32 = -X34 = Xi4 (71) 

The following operators then are symmetries of all the conserved quantities: 

= J2m nimn2m Y'^'^ = J2m "-lm"'2m 

= Em nimmm Y^'^ = J2m ^lmn2m ^^^^ 

Y'^^ = — J2m ^lmO'2m Y'^'^ = — J2m ^lrnO'2m 

They generate a fermionic realization of two copies of a Lie superalgebra: gl{l\l) x gl{l\l). The 
off-diagonal operators Y appear as the 'localized' counterparts of the bosonic operators. 



The (2, 2) — XXC model: third avatar 

There is also a 3'''^ avatar for this model where one takes A = {1,4} and B = {2, 3}. We found 
the quantities 7i^^\ C^^^ and verified that V, F and the grading are unchanged. The symmetry 
breaks down to u{l) x u(l) x u{l). However, upon localization of the fermionic generators and 
fine tuning of the twist parameters, the symmetry is enlarged to gl{l\l) x gl{V\l). Again, due 
to the boundary terms, this fermionic system is not equivalent to the previous two versions. 
Going back to bosonic variables, we find three integrable periodic-twisted boundary conditions, 
in addition to the starting bosonic periodic one. We say more about this in section (p.7|). 
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3.5 The (1,3)-XXC model 

We now consider the second and last XXC model with four states per site. This model is not 
equivalent to the above two models. It has su{3) © u(l) symmetry. As we have seen on the 
preceding two models, a particular choice of states for A and B is not restrictive, at the bosonic 
level. 

We take here A = {1, 2, 3} and B = {4}. The fermionized Hamiltonian density then reads: 



n 



mm+l — —Xi4a\^aim+in2mn2m+l — XuO.\^^iaimn2m'n2m+l 

t - - —If 

+ a2mO'2m+inimnim+l + 2:24 a2m+l'^2m«lm'^lm+l 



(73) 



^34 '3^1m^lm+1^2m+1^2ri 
^ COS 7 (1 +CmCm+l) 



2^34 0^im,+l'^lni'^2m'^2m+l 



where C = 2nin2 — 2n2 + 1. The matrices V and F, and the grading are again unchanged. The 
L-matrix reads: 



—aaim^2m 

^ '^2m'^lnT- 
Ca\^n2m 



—a aima2m 
-aaimn2m 

Ca2m,nim 



-iaa2m'nim 01^712™, 
r ■ t 

^33 —lCa2mO-lm 



ica\^a2m 



£44 



(74) 



C22 

-^33 
i244 



a n\ra'n2m -ibxu nimn2m 

-a fiimmm - ib X24 nimmm 

i a nimn2m + ibx34 nimn2m 



ianimn2m + b \x^lnimn2m - 2:^24' "'lm"-2m + i X.^lnijnn2r, 



For xi4 = X24 and the other x's unconstrained, the surviving symmetry is su(2) x u{l) x u{l). 
The generators are X^"^ and X"^^ and the four diagonal ones in (|6l|). If in additions one takes 
xi4 = X24 = —2^34, one gets a fermionic realization of the Lie superalgebra 5/(2|l), and an 
enlarged symmetry: 5/(2|l) x u{l). The complete set of generators is: X^"^ , X^^, X^^ , X^^, X^^ , 
of (|61|), yi3, of (|7|), and = n2mai^, Y^^ = Zm n2ma{^. 

There are three other avatars of this models and the analysis made for the (2, 2) model can 
be done again here. In particular, three (not four!) fermionic versions are inequivalent and we 
obtain 1 + 3 integrable periodic and periodic-twisted boundary conditions. 



3.6 The sw(4)-XXZ model 

This model, which is not an XXC one, is based on the fundamental representation of sti(4). In 
the notation of where the XXC models were further generalized, one has the (1, 1, 1, 1; 4, 4) 
system. We could fermionize the su{3) XXZ as was done for (1,2)-XXC model. We remain 
however at four states per site and study the su{4) model. 
The Hamiltonian density is given by: 



H 



mm+l 



+ "2m"2m+l 



X2'in2mn2m+1 — Xi^n2mn2m+l 



X23 n2m'n2m+l — X14 n2m?^2m+l 
2^24 ^lm^lm+1 — Xi^nimnim+l 



(75) 
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where 



F. 



(+) 



mm+l 



— X34'^a|^aim+lQ^2m+1^2m " 3^34 aima2m^2m+l 

+ COS 7 Pj^l+i + i sin 7 (P^til+i - ) 



(2) 



V-im+i't^m+i -r /i-im+i '^^m+i T /tim+l'^2m+l) (76) 
+ nimn2m{nim+in2m+l + ^lm+l"'2m+l) + "-lm^2m^lm+l?^2m+l 
= nimn2mnim+in2m+l + ^lm^2m^lm+l"'2m+l (77) 
+ ?T'l»n^2m?T'lm+1^2m+l + ^lm'^2m?T-lm+l?T'2m+l 

Let a' = sin(7 — A), 6' = sin A, c' = e*'*' sin 7 and d' = e~^^ sin 7. We find for the Lax matrix: 



mm+l 



d'aimO'2m —id'a2m'nim —'i'd'aimn2m \ 

J 



11 

-c'aim«2m ^22 -'id'a\^n2m id'al^ni 



-C aimn2r. 



£33 -id'al^ai 



c a2mni^ 



-33 

-i c'a\^a2m 



2m"'^rn 
£44 



(78) 



£11 
£22 

-^33 
£44 



a'nimn2m - b' {Xi2 nimn2m + i Xi3 nimn2m + i Xu nimn2m} 



anijnn2m + b' ^X^2^imn2m - i X23 nimn2m - i X24 nimn2mj 
i a'nimn2m + b' \^X^inimn2m - X2ini.mn2m + i X34: nimn2m\ 

24 "'\m"'2m ~r '"^34 "-lm"-2mj' 



ia'nimn2m + b' \^X^lnimn2m - a;24^"'lm?^2m + iX34 nimn2 



The grading, V and F matrices are unchanged. As the bosonic model has only diagonal symme- 
tries, these survive at the fermionic level and there are no other symmetries, whatever the values 
of the parameters x. The symmetries are the X** of (|6ll), and they generate n(l) x n(l) x u{l). 

3.7 General remarks 

Let us pause and recapitulate what we have learned so far. Several conclusions arrived at for 
the foregoing examples are of a general nature. The first concerns boundary conditions. Going 
back to bosonic variables generates many integrable systems. The non-local boundary terms are 
generically of the type El^{U^U)2 ■ ■ ■ {U^U)n~iE^ , with their hermitian conjugate. Thus, to 
each XXC system is associated a number of integrable boundary conditions equal to 1, for the 
usual periodic bosonic conditions, plus a number less than or equal to the number of versions 



the XXC model has. This generalizes to the models of [19|. The non-locality of the generalized 
Jordan- Wigner transformation appears as a strength rather than a handicap. 

Another important result concerns the symmetries. We have seen that symmetry generators 
of the bosonic model, which are local in fermionic variables, are also symmetries of the fermionic 
system with periodic fermionic boundary conditions. The generators which become non-local 
are not symmetries of the fermionic model. However we can consider their 'localized' versions, 
obtained by a simple deletion of the non-local string stemming for the fermionization. These 
odd parity operators are symmetries provided the twist parameters are chosen in a certain way. 
This also generalizes to the models of 
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All trigonometric models admit a rational limit. The rational limit of L-R is obtained by 
letting A 7A, dividing by sin 7 and taking the limit 7 — > 0. These manipulations conserve all 
the properties of the L, i?-matrices and can be applied directly to the fermionic form, without 
any other change. However the eventual extended symmetries which some bosonic models may 
enjoy breaks down for their fermionic versions. 

The general form of Vm is obvious: 



exp 



ITT 



m—1 d 



k=i j=i 



(79) 



The four phases ±1 and =bi arise from the braiding relations between Vm and the other fermionic 
operators. The particular structure of Vm is fixed by the choice of representation of C2d- This 
is true also for the grading and the matrix F. The matrix elements of the diagonal matrix F 
are found by solving a redundant system of very simple linear equations. It is in fact enough to 
solve part of this system, although we have solved it completely for the systems studied here. 

Most of these general results extend straightforwardly to the generalized Hubbard models 
considered in the following section. 



4 Hubbard models 

The Jordan- Wigner transformation turns the Hubbard model into a spin-chain of two coupled 
XX chains [Q]. In |^ ^ the construction of bosonic Hubbard models was generalized to all 
XX models. The latter are XXC models at their 'free-fermions' point: 7 = f • The fermionic 
integrable structure of the original Hubbard model was obtained by generalizing the method 



used for the XXZ model. We first briefiy recall this method |11] and then generalize it to 
fermionize generalized Hubbard models by fusing some of the models obtained in the preceding 
section. The left and right copies are labeled ] and |, respectively. 

4.1 The Hubbard model 

The Hubbard Hamiltonian is that of a spin-i electron hopping on a lattice with an on-site 
Coulomb interaction. There are four possible states per site and: 

'^2 = -J2 J2 «L+i«sm + x~^almasm+i) + U ^(2n|„ - l)(2n|„ - 1) (80) 

m s=1,l m 

Its bosonic version Q is given by 

m m 

where cr* and r* are two commuting copies of Pauli matrices. Note that this model is integrable 
without any constraint on the parameters Xs- Its L-matrix is given by 

L{X) = I{h)L^{X)CSLi{X)Iih) (82) 

where 

h = h{X) = ^Arcsinh[C/sin(2A)] and I{h) = exp(^aforfo) (83) 

The relation between h and A is the same for all models of this section. The ^/-matrix is defined 
as follows: 

Vm = V^m «> Vim (84) 



16 



where 



and fo 



Vt 



V 



and 



I'm 



Im 



^t^+l^im 



(85) 



exp 



(f Er=-/(n., 



1) j , s =t,i- The factors I{h) are diagonal and act in the 
auxiliary spaces, and are therefore not fermionized. The fermionized £-matrix reads 

Cm = I{h)C^m^£imI{h) (86) 

-^Tm — ^ -^Tm S ^ ; ~ ^ i ^ ^^'''^ 

with a grading defined by P{1) = and P(2) = 1. Both and C^m matrices are copies of 
matrix (^) where 7 = 7r/2. Note that S" is determined up to an overall factor; we take Su = 1. 
The 7^-matrix is also obtained through a similarity transformation, 

-1 5, 



n{Xi,X2)=F-'R{Xi,X2)F 



where 



F = diag(l, 1, i, i, i, i, 1, 1, -1, -1, 



-i, -1, -1) = (g) 



/I 








M 





i 














i 





Vo 








1 ; 



0h 



The corresponding grading for relation ( |30|) is given by: P(l) = -P(4) = and -P(2) = P{3) = 1. 
Two spectral parameters appear instead of just a difference because the i?-matrix of the Hubbard 
model is not additive; this is an unusual feature of all the Hubbard models considered here. 

The integrable hierarchy defined by the Hubbard Hamiltonian is known to have an 50(4) 
group symmetry (2^, 21, 22, 23, The two n(l) generators are inherited from the two 

fermionic constituents. For x-^ = = ±1, the symmetry further extends and the generators of 
the Lie algebra so(4) = su{2) x su{2) are given by: 



^ a) 

m"'lm ' 



7] 



11 



(90) 
(91) 



This extended symmetry arises after coupling of the two XX models. It exponentiates to the 
50(4) group symmetry and characterizes many physical features of the one-dimensional Hub- 
bard model (see for instance |p6|| ). 

The fermionization of the bosonic Hubbard models of is carried out in a similar manner. 
We give below the results for the fusion of some of the models studied in the preceding section. 

4.2 The (1,1) X (2,2) model 

We use as right copy the first version of the (2, 2)-XXC model. The resulting model has a local 
space of dimension eight. The Hamiltonian reads: 



n2=nl+ni + C/^(2nt^ - l)(2n^„ - l)(2n^2m - 1) 



(92) 
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where li^ corresponds to ( |20| ) and Ti^ to (|55D, both at 7 = 7r/2. The ^-matrix is a tensor 
product of V-^m of the preceding section on the left, and of 



V^N+lVlm 


\ 

' . m-1 2 














i=l j=l 

on the right. This yields for the Lax matrix C: 



where 



Cr. 



h 



I{h) C'^rn®C,imI{h) 



(93) 

(94) 
(95) 

(96) 
(97) 



I{h) = exp (^-(TfoTioj and 7^ = diag(l, 1, -1, -1) 
has been defined in the preceding section while £|m is a copy of (|57|). The grading above 

(98) 



is: 



Pt(l) = 0, Pt(2) = 1, n(l)=n(2) = 0, P^(3)=Pi(4) = l 
The 64 X 64 matrix F is given by 

F = 0-^ diag(l, i, -1, -i, i, 1, i, 1) I4 (99) 

and the grading here is: 

P(l) = P(2) = P(7) = P(8) = , P(3) = P(4) = P(5) = P(6) = 1 (100) 

This grading is used in the supertrace. 

For all x's on the right equal to each other, the symmetry of this Hubbard model is m(1)| x 
su{2)^ X su{2)i X u(l)|. Note that is is possible to take the second (or third) copy of the (2, 2) 
model. The Hamiltonian ( p2[ ) is modified accordingly. By choosing the x's on the right as 
indicated in section (|3.4[) , one restores the symmetry to: m(1)| x gl{\\\)i x g/(l|l)j. It is also 
possible to choose the model (1,3) instead of (2,2), and the appropriate choice of j;'s, with a 
resulting symmetry of: u{l)^ x gl[2\\)^ x u{l)^. The S, V, F matrices and the gradings are 
unchanged. 

4.3 The (2,2) x (2,2) model 

Again we use two copies of the first version of the (2, 2)-XXC model. The dimension of the local 
Hilbert space is sixteen. The Hamiltonian is now given by: 

n2=nl + ni + U ^(2n|i™ - l)(2n|2™ - l)(2n^™ - l)(2n|2m - 1) (101) 

m 

where ^2 and corresponds to two copies of (^) at 7 = 7r/2. The numerical coupling matrix 
is given by: I{h) = exp(|C|oC'|o) where C = 7^ = diag(l, 1, —1, —1). The ^-matrices are: 



exp(f EI^I'E. 



j=i{nsji 



1) , 



4,i 



(102) 
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The Csm matrices are copies of matrix ( p7| ) at 7 = tt/2. The grading in (86-8^) is given by: 



P(l) = P(2) = 0, P(3) = P(4) = 1, and S = diag(l, for the 'up' copy. For (|3g) and the 
supertrace the grading is: 

P(l) = P(2) = P(5) = P(6) = P(ll) = P(12) = P(15) = P(16) = 
P(3) = P(4) = P{7) = P(8) = P(9) = P(10) = P(13) = P(14) = 1 

The 256 x 256 matrix F is given by 

F = o-'' (g) I2 O diag(l, -1, i, i, -1, 1, -i, -i, i, -i, 1, 1, i, -i, 1, 1) I4 (104) 

For all x's, on both left and right independently, equal to each other, the symmetry of this 
Hubbard model is su{2)i x sm(2)| x u{1)i x su{2)i x su{2)i x n(l)|. As indicated at the end of 
section ( [4.2D , it is possible to fuse any two copies of the other versions of (2,2) , or (1,3). The 
resulting symmetry changes accordingly, with the appropriate choices of x's. The gradings and 
V, S, F are unchanged. It is possible to couple also (1,2) models, or in general any two XX 
models in their fermionic form. A fermionic Hubbard Hamiltonian is simply: 

n2 = nl + ni + uY,clci (105) 



4.4 Symmetry and Integrability 

There are many models of itinerant electrons in the literature which qualify as such for the 
name Hubbard model. However, most of these models have an additive i?-matrix, or do not 
share the same algebraic structure of the original Hubbard model (the Bariev model). Here we 
chose to use the term Hubbard models for a class of models which share the same integrable and 
algebraic structure as the original model along with a similar structure for the non-additive R- 
matrix. These bosonic multistates Hubbard models were constructed and studied in |6|, |7| . (Lax 
pairs were derived in ||2^.) They correspond to arbitrary up and down copies of the generalized 
XX models, with and on-site coupling. This structure, which is the natural generalization of 
the Hubbard model one, has been fermionized here, and we obtained the first multispecies 
generalizations written in fermionic form. 

Recall that (see section ( |4.1] )) the Hubbard model has, after coupling of its two components, 
an extended symmetry not shared by its two independent components. This symmetry also 
exponentiates to a group symmetry. It is therefore natural to ask whether the new models 
we obtained here have such a property. We have looked for additional symmetries of the new 
Hubbard models, and found it is unlikely that they exist. This negative result can be explained 
by a conflict between higher symmetries and the strictures of integrability. The form of the 
fermionized XXC models shows that they are not constructed out of group invariants and their 
symmetries do not exponentiate. Coupling does not solve this problem. There is however a finite 
number of discrete symmetries which appear; they correspond to fermion species interchanges. 
In fact one could easily build models of itinerant electrons with larger symmetries. However 
integrability breaks down the symmetry group manifold to a discrete set of points, some cor- 
responding to the symmetries of the XXC blocks and some appearing only for the Hubbard 
models. 



5 Non integrable models 

The fermionized Hamiltonians seen so far were obtained directly from their bosonic counterparts 
without use of any integrability criterion. The generalized fermionization procedure is indepen- 
dent of any integrable structure. We now fermionize two non integrable models of considerable 
physical interest. 
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5.1 The spin-1 XXZ Heisenberg chain 



The bosonic Hamiltonian is given by: 



(106) 



where the S** form a spin one representation of su(2), at every site. For arbitrary the diagonal 
operator 5^ = Yl,m^rn commutes with For A^ = 1, 5^'^^ = J2m^m^ ™ addition commute 
with the Hamiltonian. 

To fermionize this model one has four choices of embeddings. We chose the same one as in 
section ( |3.2| ). All states with double occupancy on at least one site are therefore excluded. The 
fermionization then yields: 

(•^Im+l^lni + ^Im'^lm+l) ^2m^2m+l 
„t J _ , J J 



n 



mm+l 



(107) 



+ 
+ 
+ 
+ 



lm+l'^2m'^lm-'^2m+l + aim'^2m+l'^lm+l'^2m 

1) ^2m 



(o2m+l'^l"T''^lm,+l + '^lm+1 '^lm'^2m 



VTri 



-1 ('^lm+l"lm'^2?Tt + «2m"lm«lm+l) "-2771+1 



(2n2m - l)(2n2m+i - 1) 



where 

7ro = l , vr^^i = (2nii - l)(2n2i - 1) (2ni„_i - l)(2n2^_i - 1) (108) 

Note that, contrary to all the integrable models we have studied so far, there are non-local 
contributions associated with odd-parity combinations in the Hamiltonian density. It is also 
easily seen that its action is stable on states without double occupancy. Let p be the projection 
operator (pO|) defined in section ( p.2[ ). The action of 7^2 = Z^m^mm+i on the S-'^-dimensional 
subspace is equivalent to the action of p7i2P on the 4^-dimensional space. 

We may remove the vTm-i's, leaving their factors in, wherever they appear in (107), and 
obtain a local density which combines bosonic and fermionic pieces. The local spin-one fermionic 
version the XXZ Heisenberg model may then be defined this way. The diagonal operator = 
^^nim(2n2m — 1) still commutes with the Hamiltonian with periodic boundary conditions. 
However as seen for the spin-^ Heisenberg chain, the extended symmetry at A^ = 1 does not 
survive the fermionization. 



5.2 The spin- 1 XXZ Heisenberg chain 

The bosonic Hamiltonian is given by (|106| ) where now the 5* span a spin | representation of 
su{2). Here the Hilbert space of the chain is 4^-dimensional. We find for the fermionized 
Hamiltonian density: 

_ -'^ t [q t I Q t t I Q 



n 



+ 3 a2ma2m+l + 4 "'2m?T-2m+l 



(109) 



+ 



1 



3 02m+l'^2m + 3 a2m+l'^2m "I" 3 a2m+l02r 



+ 3 a2m+ia2m + 4 n2mn2m+l 



+ \/3 
A 



«2m"'2m+l — 02r)i'^2m+l — a2m+1^2m — a2m+l"'2m 
«2m"'2m+l ~ a2m^2m+l + a2„^in2m + a2m+1^2m 



+ 



{2nim - l)(4n2m - l)(2nim+i - l)(4n2m+i - 1) 
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Remarks similar to those made for the spin-1 case hold here. In particular a local density may 
be defined by the deletion of TTm-i- The diagonal commuting operator is given by 

= ? E(2^1- - l)(4n2m - 1) (110) 

6 Conclusion 

We have introduced a generalized Jordan- Wigner transformation which allows the fermioniza- 
tion/bosonization of integrable and non-integrable spin chains. Several variants of this transfor- 
mation were described. We have shown on various examples how the mapping works. We also 
obtained the first fermionic versions of generalized Hubbard models. 

The study of particular models allowed to infer some general results. In particular we showed 
that the non-conservation of boundary conditions, far from being a 'shortcoming', is in fact an 
advantage and generates many new integrable boundary conditions of the periodic-twisted type. 
Another issue is the mapping of certain local operators to non-local ones. However we have 
shown that, the non-local fermionic operators can be made local and the symmetry of the 
model turned into supersymmetry, without any loss of integrability. Fermionic realizations of 
Lie algebras and superalgebras appeared naturally in this context. 

Periodic boundary conditions were considered here. But the fermionization method can 
be directly applied to models with other types of boundary conditions, within or without the 
framework of integrability. 

Fermionization and bosonization of one-dimensional systems are an expression of a local 
equivalence between a bosonic language and a fermionic one. This is clearly seen at the Hamil- 
tonian level, for integrable and non integrable models. The fermionization of the L-matrix is a 
local procedure. The important point is to notice that a periodic boundary condition breaks this 
equivalence, while an open one, where -ffAri is dropped, does not. (Integrability may be lost but 
this is a secondary point.) The issue of boundary conditions is intrinsic to any fermionization 
scheme. 

The present work was near completion when appeared. The authors of |^] proposed a 
fermionization scheme which applies only to integrable models. Their main example corresponds 
to the one in section (|3.2| ), at 7 = tt/2. The Hamiltonians are given by the same expressions while 
^-matrices have the same structure. The phases zizi do not appear, and this is general, in the 
method of ||8|. It would be interesting and instructive to find, for integrable systems, the exact 
relation between the two methods. However let us stress that boundary issues invariably arise 
for any fermionization method, and only the foregoing approach tackles this issue. The method 
of Gohmann and Murakami and the generalized Jordan- Wigner transformation are therefore 
complement ary. 

A given one-dimensional Hamiltonian density can be transposed on a higher-dimensional 
lattice. While in one-dimension bosons and fermions are closely related, the relationship breaks 
down in higher dimensions. This is one important motivation for finding the fermionic expression 
of a bosonic model in one-dimension, before going to higher dimensions. The physical properties 



of a model may change drastically with dimension. As noted in section (4^), the study of one- 
dimensional models would shed some light on the interplay between integrability and symmetries. 
But it is only on a two-dimensional square lattice with nearest-neighbor interaction that the 
Hubbard model exhibits d-wave superconductivity. It would be worthwhile to study the new 
Hubbard models in one and two dimensions. 

Acknowledgement: Z.M. would like to thank N.E. Bickers for sharing some of his insight 
on the Hubbard model. 
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